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! Abstract 

The gluon fields screening in the stochastic vacuum of gluody- 
namics is studied. The effective action is derived for the instanton 
interacting with nonperturbative fields. Quantum nonperturbative ef- 
fects are shown to affect greatly the shape of instanton. The power 
asymptotics x~ 2 of the classical "instanton's profile function" at large 
e-j ! distances is replaced due to yhese effects by Airy function asymptotics 

x~ l l i exp{— const ■ x 3 ^ 2 ). 

1 Introduction 

In the last years a systematic description of nonperturbative effects in QCD 
has been provided in terms of the gluon nonlocal gauge-invariant correlators 
[1]. This Vacuum Correlators Method turned out to be very successful in phe- 
nomenological description of important QCD phenomena [2]. By now, there 
exist numerical results from lattice simulations concerning the fundamental 
field strength correlators for pure gauge theory with gauge group SU(2) [3] 
and SU(3) [4] over physical distances ranging up to 0(1) fm. 

On the other hand during the last 15 years we have witnessed active 
development of nonperturbative QCD picture based on the instanton liquid 
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model [5,6,7]. It is therefore of great interest to compare the lattice calcula- 
tions of field strength correlators with the computation of the same quanti- 
ties in instanton-anti-instanton vacuum model. The analytical calculations 
of field strength correlators for the instanton-anti-instanton vacuum in the 
dilute gas approximation exhibit the power decrease of the bilocal correlator 
at large distances. On the other hand, the results of lattice simulations are 
consistent with exponential decrease of this correlator. Thus we see that the 
computation of the field strength correlator requires the comprehensive anal- 
ysis of instanton field character in nonperturbative vacuum. It is necessary 
to note that already in first publications on the instanton liquid model the 
Debye screening was revealed with a typical mass of the order of 350 MeV 
[6]. A new mechanism of instanton scale stabilization essentially different 
from that considered in [6] was proposed in [8,9]. It is due to the effective 
"freezing" of strong coupling constant in the nonperturbative vacuum. As a 
result the Debye screening mass becomes smaller. 

In the present paper are obtained the effective action of topologically 
nontrivial fluctuations in the nonperturbative vacuum making use of the 
Vacuum Correlators Method. The new type is revealed of instanton screening 
different from previously considered [6,10]. 

This paper is organized as follows. In Section 2 the effective action is 
derived for the instanton-like object in the nonperturbative vacuum. In 
Section 3 arguments based on scale analysis are presented to show that this 
effective action which breaks scale invariance possesses nontrivial solutions 
stable in size. In Section 4 are investigated the solutions of the equations of 
motion in asymptotic regimes and the instanton field screening is studied. 

2 Effective action 

The Euclidean action of gluodynamics is written as 

S = ± j #xtrJ%{x\ (1) 
where we use the Hermitian matrix form for gauge fields is used, 

•Ffiv = d^Av d u J\.^ ^[.Aju, -4-1/] , -A-fi = QA^t > trt t = —5 (2) 
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To obtain effective action for instanton in the nonperturbative vacuum we 
represent the full field A^ as 

A li = A li + (3) 

where is the external field (instanton-like object with unit topological 
charge) and are nonperturbative vacuum fields. 

Under gauge transformations (U is the transformation matrix), the fields 
change as 

A„{x) -> U + (x)A„(x)U(x) 

B»^U + {x){B^x)-id»)U(x) (4) 

The general expression for the effective action of an instanton in the 
nonperturbative vacuum has the form 

Z = e S eff [A] =< e S[A + B ]+ S[B] >b> (5) 

where 

< 6(B) > B = J dp[B]6(B) (6) 

and dfi[B] is the measure of integration with respect to nonperturbative 
fields. The explicit form of this measure will not be needed in our analysis. 
The full field strength can be rewritten in the form 

F, V [A + B] = F^[A] + G^[B] - i[A„, B v ] - i[B„, A,] (7) 

Let us now consider a large-size instanton. The squared strength of the 
instanton field at the center is given by 

W(* = *o)) 2 = 192/p 4 (8) 

For p > po = (192/ < G 2 >)^ 4 ~ Ifm, the instanton field is weak relative 

to the characteristic field strengths in the vacuum gluon condensate 

< G 2 >=< (gG^) 2 >~ 0.5Gev 4 [11] and can be treated as a perturbation. 

Therefore, expanding the effective action (5) up to the second-order terms 
in A^ we have 

S eff [A] = -InZ ~ S[A]+ < S dia [A, B]> B + < S para [A, B] > B , (9) 

where 

S dm [A,B] = ~ J d^xtr{[A^B v \ + [B^A V ]} 2 , (10) 
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S P ara[A, B) = J #xtfytrF^(x)G^(x)frF aX (y)G aX (y), (11) 

and the conditions < B^ >= 0, < G^ v >= are taken into account. 

Hence, there exist two contributions [9]. First, there is the motion of 
"charged particles" in the external gluonic field F^. This motion leads to 
screening similar to orbital Landau diamagnetism. And second, there is the 
direct interaction of F^ u with the spin of the field B^. In the second order 
in the external field, this interaction yields an antiscreening term, which is 
analogous to that describing the Pauli paramagnetic effect. 

A similar decomposition into diamagnetic and paramagnetic terms was 
suggested by Polyakov [12] for perturbative gluodynamics, and led to the 
simple and elegant explanation of antiscreening with correct coefficirent for 
/^-function. 

Let us consider paramagnetic component of S e ff. The general form of 
the instanton field is 

A^x) = Xq)VlA^(x - x ))fi + $(x , x), (12) 

where 

rx 

$(x,x ) = P exp{i / Bpdzn] (13) 

Jxo 

is the operator of parallel transport, Q is the matrix of global rotations in 
the color space. A^(x — Xq) in the equation (12). 

4 = ^% v %f{A (14) 

is the instanton field in some gauge. The / is so-called "instanton's profile 
function" . Classical function / has the standard form in the singular gauge, 

fd = TTF77) (15) 

Prior to averaging over the fields B^, we write down the paramagnetic term 
in S e ff in the following formQ 

S para [A,B] = -— J d A xd A y J dfltrflF tlu (x)fl + G fll/ (x,x ) 

xtrnF aX {y)tt + G aX {y,x ), (16) 



1 In this section the notation = for brevity is used. 
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where 

G^x, x ) = $(ar , x)G llv {x)^{x, x Q ) (17) 

and dVt is the Haar measure on the SU(N C ). In the tensor notation, the 
integrand in (16) is written as 

^a l p 1 (F^)p 1 a 2 ^a2P2(G^)p 2ai ^ a3 p 3 (F a x)/3 3a4 il^p i (G (7 x)p i a3 (18) 

Taking into account the relation 



/ 



dQQa-^fa fi^ 2/32 ti a3 p 3 ^a 4 /3 4 



and integrating (16) over the Haar measure, one obtains 

s para [A,B] = ~ 2 (7v 2 - i) g 4 / d4xd4 y trF ^( x ) F ^(y) trG ^( x ^ x o) G ^(y' x o) 

9 (20) 
Let us perform averaging over the nonperturbative field and remind that 
the nonlocal gauge-invariant correlation function < trG^ u (x, xo)G a \(y, xo) > 
can be represented as [5] 

< G 2 > 

< trG^(x, x )G aX (y, x ) >= — — — {(5^5 uX - 6 lt x6 va -)D(x - y) + 0{D 1 )} 

(21) 

In this way, the paramagnetic term in S e ff is reduced to the form 

S para [A] =< S para [A,B] > B = - ^ N T_\ )gA J d 4 xd%F^(x)D(x-y)F^(y) 

(22) 

Let us consider the diamagnetic part 

S dm [A,B] = --J d A xtr{[A^B v f + [A,,, B U ][B„ A,]}. (23) 
9 J 

Again we integrate over the Haar measure first. Thus, the Sd ia is given by 
the following expression: 

S dia [A, B] = ^^%r / d'xtrAl{x)trBl{x) (24) 



To obtain gauge-invariant expressions for vacuum correlator, we use the 
Fock-Schwinger gauge x^B^ = . By virtue of the fact that the t'Hooft 
symbols f\ a are antisymmetric, the instanton field satisfies this gauge condi- 
tion automatically. In this gauge, the nonperturbative field _B M is related to 
the field strength G^ u by the well-known equation 

Bfi{x) — x v I adaG^(ax) (25) 

J 

Using the relation (21) we arrive at 

S dia [A] =< S dia [A,B] > B = ± J d 4 xJ:(x)(A;(x)) 2 , (26) 

where 

= -^ N2 C _ 1 < G 2 > x 2 ^ ada ^ (3d(3D(\a - (5\x) (27) 

Bringing all the results together, we find that S e ff for instanton-like object 
in the nonperturbative vacuum is given by 

S eff [A] = S[A] + S para [A] + S dta [A] 

= 4^ / * X *VW< X ~ V)W + d 4 xZ(x)(A°(x)) 2 , (28) 

where 

e(x - y) = ^(x - y) - U(x - y), U(x - y) 
Equation of motion for A^ is 

Vf(A(x))Jd 4 ye(x-y)Fi l/ (y) 

and 

V)f = 5 ab d^ - if abc A c ^ 
is the covariant derivative on the gauge group SU(N C ). 



= mf^w D(x - y) (29) 

- E(x)A a „(x) = 0, (30) 
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3 Scale tranfor mat ions 



In the previous section we have derived S e ff for weak external field (instanton- 
like object of large size) in the nonperturbative gluon vacuum. It is known 
that scale invariance in gluodynamics is broken due to the anomaly of the 
trace of energy-momentum tensor < 9^ >~ — < G 2 >. Topology allows 
the existence of topologically nontrivial solutions. However their stability 
in size is not evident from scale arguments. For example in theories which 
include Higgs fields (Weinberg-Salam or Georgi-Glashow models) instanton 
solutions are absent at the classical level. Instanton tends to contract in size 
thus lowering the classical action and finally becoming a point singularity 
with p = 0. 

If however one considers the effective action S e ff with quantum correc- 
tions over the instanton background field, one can perform a scale transfor- 
mation x^ — > Xfj,/X, A^(x) — > X^A^x/X) and then consider the derivative 
of S e ff with respect of the logarithm of the transformation parameter. This 
yields 



where 0^^(x) is the trace of energy-momentum tensor for a given field config- 
uration. The trace 0^(x) contains a positively defined contribution which is 
due to the mass term of the classical action. This mass term breaks scale in- 
variance and leads to the collapse of the classical instanton. In the quantum 
theory 6^ acquires a negative contribution due to quantum anomaly. It is 
the presence of this anomaly which eventually stabilizes the instanton. One 
may argue that the size of the instanton is determined from the condition of 
the total pressure to be zero [13] 



Let us consider the scale transformation of the effective action (28) and 
prove the existence of the stable in size instanton in nonperturbative scale 
noninvariant gluodynamics. 

Numerical lattice calculations revealed that, to a high accuracy, the func- 
tion D can be approximated as [3,4] 




(31) 




D(x-y)= e-ri x - v \ [i = 1/Tg ss IGeV. 



(32) 
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Let write (22) in the Fourier form 



= -A/ 704^(?)^(-9)n(g 



4g 2 i (2tt 
Taking into account (29) one obtains 
<G 2 > 



<G 2 > 



7T 2 // 



12(iV2 - l)g» 

Perform now a scale transformation 



(33) 



(34) 



(iV c 2 -l)^ 2 (/i 2 + g 2 ) 5 /2 
Ag M , - A 2 F; i/ (Ag). Then 



95, 



para 



dlnX 



A=l 



1 



/ J^Lp a (q)F a ( 



(2vr 



an(g 2 ) < Q 

dlnq 2 ~ 



(35) 



Consider in the same way. Performing a transformation 



x,. 



dSdia 



dlnX 



A M^ir/A) one arrives at 
^ / d 4 x(^(a;)) 2 (2S(x) + 



dlnx 



(36) 



Making use of the relations (27) and (32) and of the conclusions on 
D(\a — (5\x) from Ref. [9] one can deduce that 2S(x) + dT,(x)/dlnx > 
Therefore 

dSdia 



dlnX 



> 



(37) 



A=l 



It means that diamagnetic interaction of the instanton with nonperturbative 
fields leads to the collapse of the instanton. On the other hand the paramag- 
netic interaction effectively results in instanton "swelling" in its scale. Hence 
a stable in size instanton can exist in the nonperturbative gluonic vacuum de- 
scribed by nonlocal gauge-invariant correlator < trG flu (x,x )G a \(y,x ) > . 



4 Screening 

The instanton profile function in the nonperturbative vacuum is determined 
by Eq. (30) for the ansatz (14). At small distances from the center of the 
instanton x — > 0, (q — > oo) equations (27), (29) and (34) lead to | 

E(z - 0) -> < G 2 > x 2 

V ; 16 iV 2 - 1 



2 In this section, we use the notation x — \ x 1 
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e(q -> oo) -> 1 (38) 

So at small distances one retrieves the classical equation for the profile func- 
tion 

% + 2 fx ~ l (2f ~ 3/2 + /} = °' /(0) = 1j /(00) = ° (39) 
and, correspondingly, f(x — > 0) — > / c ;(x) (15). 
At large distances x — > oo, (g — > 0) 

E(x - oo) ^ < G 2 > T g x(l + O(^)) 

r 2 



e ( 5 _> 0) ~ 1 - 7-^-^ < G 2 > T 9 4 (40) 



7T 

{Nf — l)g 2 

Thus, we have the following asymptotic equation for the profile function 

0- 7 x/ = O, /(oo) = 0, (41) 

where 

7 4^< G! > T '< 1 -(7vr ! IF <G2>T » rl (42) 

Hence the function / is asymptotically damped at large distances like the 
Airy function 

f(x) ~ Ai(^ 3 x) - I./l 7 -Vi2 a: -i/4 ea; p(_2 1/2a _ 3/2 a; — > oo (43) 

2 V 7T 3 

One notice that from (43) we have screening of the instanton profile function 
with effective mass m sc = (2/3) 2 / 3 7 1//3 xs 7 1 / 3 . Choosing the standard values 
of the parameters in QCD N c = 3,T g = lGeV' 1 , < G 2 >= 0, 5Ge\/ 4 , a s = 
g 2 /47r = 0,3 one obtains the screening mass m sc (QCD) ~ 250MeV. On 
the other hand, as it is well-known, in gluodynamics the value of the vac- 
uum condensate is < G 2 >gd— 4 < G 2 > [14]. Therefore , one obtains the 
screening mass m sc (GD) ~ 520MeV. 
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5 Conclusion 



One concludes that the effect of screening is able to change drasticully the 
interaction in the instanton liquid, where the average distances between pseu- 
doparticles (~ (GOOMeV)^ 1 ) are of the same order of magnitude as inverse 
screening mass (~ (500MeK) _1 ). The numerical calculations for the screen- 
ing effects in the instanton vacuum and comparison with the lattice calcula- 
tions are in progress. 
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